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Abstract 



The accuracy of the conventional treatment of quark mass thresholds in the QCD running coupling 
pH ■ constant a^, based on the step approximation to the /3-function, is investigated. The errors of 

a , extrapolating as(/i) from low energies to ji = Mz, implied by this approximation, are shown to be 

of the same magnitude as the typical next-to-next-to-leading order contributions to ag and tend 
to increase the resulting Og ■ The importance of the proper choice of matching points is emphasized. 



Q>^ I In this note I shall discuss the importance of the proper treatment of quark mass effects on the QCD 

OO ■ running coupling constant . Although the quark mass effects are not large, steady improvement 

in the precision of experimental data, coming in particular from a new generation of experiments 
I at CERN and Fermilab, combined with significant progress in higher order QCD calculations, have 

' led to a renewed interest in quantitative aspects of these effects fl], . There are essentially two 

■ reasons for it. 

1 ^ ' The first concerns the exploitation of the next-to-next-to-leading-order (NNLO) QCD calcu- 

I ' lation that have recently become available for a number of physical quantities 1^, |^, ^, ^ and which 

exist basically for massless quarks only. The NNLO corrections are tiny effects and to include them 
r~| , makes sense only if they are large compared to errors resulting from the approximate treatment of 

the quark mass thresholds in massless QCD. To quantify the importance of the NNLO correction to 
the couplant as, consider the difference between the values of as(MS, Mz) in the NLO and NNLO 
approximations, taking, as an example, A^^i- = 0.2 GeV. At the NLO the couplant a = cks/tt is 



X 

5t , given as a solution to the equation 



where for Uf massless quarks 



61n-!^ = - + cln (1) 

A(5) a l + ca ^ ^ 



33-2n/_ _ 153- 19n/ 
6 ' ^~ 33 -4n/ 



are the first two, universal, coefficients of the QCD /3-function appearing on the r.h.s. of the 
definition equation for the couplant a(/i) 



da{iJ,) 
dln^ 



I3{a) = -ba'^ifi) (l + caifj.) + C2a^{fi) + •••). (3) 



For ij = Mz,nf = 5 and = 0.2 GeV we get a^^'^{MS,Mz) = 0.03742. At the NNLO and for 
C2 > c^/4 we have instead of (||) 

Mn ^ = i + cln '\ + /(a, c^) (4) 

A a Vl + ca + C2a^ 
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2c2 — (? ( 2c2a + c c 



f{a, C2) = l^arctan arctan -J ; d = \J 4c2 — (? (5) 

where C2(MS,nj = 5) = 1.475. Combining % and (|) yields a^^^o^j^Js^ ^ 0.03665. The 
relative difference between these two approximations 



^NLO _ ^NNLO 



a 



NLO 



0.02 (6) 



thus amounts to about 2%. The NNLO correction should therefore be included only if the neglected 
effects can be expected to be smaller than this number. However, as we shall see, quark mass 
threshold effects can in some circumstances be just of this magnitude! In general the magnitude 
of higher order corrections to depends on the renormalization scheme employed and so does 
also the estimate (P). However, if defined as the relative difference between the NLO and NNLO 
approximations this dependence is weak. 

The second reason is related to the problem of comparing the values of a<j, determined from 
different quantities characterized by vastly different momentum scales. As recently emphasized in 
an extensive review on determinations Q, there is small, but nonnegligible discrepancy between 
the value of as(MS,M^) obtained by extrapolation from some of the low energy quantities and 
as(MS,M^) determined directly at the scale Mz at LEP, the latter giving the value higher by 
about 5-10%. Simultaneously, it has been noted in that there is an exception to this behaviour 
in the case of the ratio Rr which, when extrapolated from ra-j- to Mz-, gives values of close to 
that measured directly at LEP (0.120 it 0.005 |9|) as well. The physical relevance of the discrepancy 
between the low energy extrapolations and direct measurements of at LEP has very recently 
been emphasized by Schifman in ||l^. In particular, he argues that the extrapolations of asivn^^ 
to as{Mz^ are unreliable due to problems with the control of power corrections. As will be shown 
below, a part of this overestimate of the extrapolated value of as(M^) may in fact be simply due 
to the approximate treatment of the c and h quark thresholds. 

In the rest of this note I shall analyze the quantitative consequences of the exact treatment 
of quarks mass thresholds at the LO and formulate the conventional matching procedure [0] for 
massless quarks in such a way that its results are so close to the exact ones that the available NNLO 
calculation can be consistently included. I shall describe in detail the approximation in which the 
"light" quarks li, d and s are considered massless, while the c, h and t quarks remain massive. 

As complete multiloop calculations with massive quarks are very complicated and are available 
only at the leading order, all higher order phenomenological analyses use the calculations with a 
fixed effective number of massless quarks, depending on the characteristic scale of the quantity. 
In order to relate two regions of different effective numbers of massless quarks the approximate 
matching procedure developed in [11| is commonly used. It should be emphasized that this pro- 



cedure concerns only those mass effects that can be absorbed in the renormalized couplant. At 
higher orders there are, however, mass effects that remain in the expansion coefficients even after 
the effects of heavy quarks have been absorbed in a suitably defined running couplant. 

In QCD with (at least some) massive quarks the renormalization group equation for the couplant 
a(;u) looks formally as in massless QCD. The only, but important, difference concerns the two lowest 
order /^-function coefficients, 6, 6c, which are no longer unique as in massless QCD, but may depend 
on the scale /i. While in the class of MS~like renormalization conventions h = 11/2 — n//3 as in 
massless QCD, in MOM-like ones it becomes a nontrivial function of the scale ji |jl2| 

h{ii/mi) = -]-^hi{xi), Xi = — (7) 
where the sum runs over all the quarks considered, are the corresponding renormalized quark 
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Figure 1: a) 6(/i/mj) together with three step approximations, corresponding to matching at the 
points /ij = Kirii'ji = c,b,t with k = 1, 1.58, 2.24; b) the ratio Ra for the same three approximations, 
made to coincide at //q = 1 GeV. 



masses ^ and the threshold function h{x) is given as |12 



h(x)=6x^ dz V— J— = 1- — + In-p = ^. 8) 







The last - approximate - equality is a very accurate approximation to the exact form of h{x) in the 



whole range x S (0, oo) |12|. This allows a simple treatment of the quark mass thresholds at the 
LO. There is no analogous calculation of the next /3-function coefficient, c, for massive quarks and 
thus no possibility to repeat the present analysis at the NLO. This is one of the reasons why most 
of the phenomenological analyses use the so called "step" approximation in which at any value of 
fj, one works with a finite effective number of massless quarks, which changes discontinuously at 
some matching points //j. Consequently b{nf) becomes effectively a function of /i, discontinuous at 
these matching points, as shown in the Fig. la. 

The matching points are assumed proportional to the masses of the corresponding quarks, 
fii = Kirii. In principle different quark threshold could be associated with different k, but for 
simplicity I take them equal. The free parameter k, allowing for the variation of the proportionality 
factor, turns out to be quite important for the accuracy of the step approximation. At the LO the 
matching procedure then consists in the following relations at the matching points /Uj (the number 
in the superscript defines the corresponding effective number of massless quarks) 

/ .(3)\1/3M4) 

a]SA-mJK^^y) = a^O'^Km./A(^)) ^ = A(3) A_ (9) 



1/36(5) 

a]SA^m,/K^^) = a^o.5(^^,/A(5)) ^ A(5) = A(^) ( ^ ) (10) 



KTUb 



^For the purpose of this discussion quark masses are regarded as constants, altough they too get renormaUzed and 
thus "run" as well. 
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/ .(5)\l/3''(6) 

a^O,5(^^^/A(5)) ^ a^O'6(^^^/A(6))^A(6)=A(5) -— (11) 

Note that each of the intervals of fixed n-j is associated with a different value of the A-parameter, 
A("/). The resulting dependence a(/x/mi) on /x is thus continuous at each of the matching points, 
but its derivatives at these points are discontinuous, reflecting the discontinuity of the step ap- 
proximations to 6(/u/rn-j). This procedure can be easily extended to any finite order. Let me point 
out that the more sophisticated procedure for matching the couplants corresponding to different 



effective rij developed in ||l^ coincides to the LO with (^)-(ll). 

To estimate the errors in resulting from the above defined approximate treatment of quark 
thresholds we merely need solve the LO equation with exact explicit mass dependence as given in 

m 



da{fi) 



-a' 



(f-^m. (12) 



dln/i 

For our purposes the approximation h{x) = x^/(5 + x^) is entirely adequate and yields 

a(/x) = ; (13) 



11 3\ , /i 1 V- 1 V/i^ + 5m2 
in — 7— > m ■ 

A (3) q /-^ 



2 3 J A(3) 3 . ^ , //.nu2 



i=cAt V^(A{3))^ + 5m2 



where the fraction | comes from the sum over the three massless quarks u, d, s and A^^^ is the 
corresponding A-parameter appropriate to 3 massless quarks. For the heavy quarks c, b and t I 
take in the following nic = 1-5 GeV, = 5 GeV, rrit = 170 GeV. The distinction between the 
"light" and "heavy" quarks is given by the relative magnitude of mi and A, the latter being defined 
by the condition 5m? ^ A. For the above values of mc,mi,,'mt this condition is very well satisfied. 
Consequently, for /i <C rrii, i = c,b,t (|l^) approaches smoothly a^'^ for n/=3, while for ^ ^ rrii and 
neglecting A^^^ with respect to 5m?, it goes to 

«(/^) = T-^ 7= 7= — ^ = ^—n (14) 

h(f^-)]n ^ I 1 . / V V5mfc V5mf \ &(6) In ■ 



A(3) A(3) J ^^'HV5) 
where the parameter A^^\k) depends in general on k and 



(15) 



coincides with A^^^ defined via the subsequent application of the matching relations (P)~(|Tl|) for 
K = y/b = 2.24. Although from the point of view of the matching procedure k is not exactly fixed, 
the value k = ^/b will be shown to be in some sense the best choice. The relation between A^^^ and 
A^^) depends nontrivially on k. 

In Fig. la b{^) is plotted as a function of ^ for the above mentioned masses of c, h and t quarks, 
together with its step approximations corresponding to three different values of «; = l,^/b, ^b/2. 
There is hardly any sign of the steplike behaviour of the function h{x) in the region of the c and 
b quark thresholds and only a very unpronounced indication of the plateau between the b and t 
quark thresholds. The steplike approximations are poor representations of the exact h{x) primarily 
due to the rather slow approach of h{x) to unity as x — > oo. However, there is a marked difference 
between the three approximations. While the step approximation with the conventional choice 
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K = 1 understimates the true h(x) in the whole interval displayed and would do so even when some 
smoothing were applied, k = \/5 gives clearly much better approximation as the corresponding 
curve is intersected by the exact h{x) at about the middle of each step. 
In Fig. lb the fi dependence of the ratio 



between the above exact solution (14) and the approximate expressions for the above mentioned 
values of k is plotted assuming A(3~ 200 MeV 0. As we basicaly want to compare the results of 
different extrapolations starting from the same initial /xq, A^^^ used in the approximate solutions was 
rescaled by the factor 1.004 with respect to A^'^^ in (p!3|), thereby guaranteeing that all expressions 
coincide at //q = 1 GeV. Any deviation from unity in Fig. lb is then entirely the effect of an 
approximate treatment of the heavy quark thresholds. The Fig.l contains several simple messages. 
The approximate solutions based on the matching procedure defined in (^)-(ll) are in general 



much better immediately below the matching point than above it and is worst at about Smmatch- 
This reflects the fact that the function h{x) vanishes fast (like x^) at zero but approaches unity only 
very slowly. Moreover, in the Mz range the effect of the c quark threshold is essentially the same 
as that of the b quark and both are much more important than that of the top quark, although 
Mz/rric ^ 60, Mz/mt « 18, while Mz/mt « 1/2! 

The effect of varying k is quite important, in particular with respect to the c and b quark 
thresholds. In general k > 1 improves the approximation above the matching point, but worsens it 
immediately below it. The choice k = \/5, suggested by the asymptotic behaviour of (|l^ ) is clearly 
superior practically in the whole displayed interval /i G (1, 10^) GeV and leads to an excellent 
(on the level of 0.1%) agreement with the exact solution in this interval. On the contrary the 
conventional choice k = 1 leads to much larger deviation from the exact result, which exceeds 2% 
in most of this region. This discrepancy is thus of the same magnitude as the effects of NNLO 
corrections to the couplant itself. It turns out that the effect of an exact treatment of the quark mass 
thresholds is as important as that of the NNLO correction to the /3-function and must therefore 
be taken into account, whenever the latter is considered. 

In Fig. 1 the exact as well as approximate expressions for the couplant are made to concide at 
^0 = 1 GeV, in order to stay below the c quark threshold. In actual extrapolations one usually 
starts at somewhat higher scales. For the r-lepton decay rate as well as the GLS sume rule, the 
corresponding initial fiQ is about 2 GeV, for the various T characteristics /^o ~ 10 GeV. In Fig. 2 
the ratio (|l6|) is therefore plotted for these two starting scales and the indicated values of k. We 
see that 

• For /io = 2 GeV the situation is essentially the same as for = 1 GeV, i.e. extrapolated 
a^'^{Mz) is about 2% higher when the step approximation with k = 1 is used. The fact that 
conventional extrapolations of asimr) to as{Mz) use this step approximation with k = 1 



could thus cause part of the overestimate of as{Mz) discussed in [|lO|. On the other hand the 
choice K = \/5 again leads to an excellent agreement with the exact extrapolation above 10 
GeV and very satisfactory even down to 1 GeV. 

For fiQ = 10 GeV the choice k = \/5 is only marginally better that k = 1 when extrapolating 
to /u = Mz, but now both values yield a^'~'{Mz) only about 0.5% from the exact LO result. 
This observation is consistent with the fact that as{Mz) extrapolated from T characteristics 
is lower than that from Rr and consistent with those based on deep inelastic scattering. For 
this choice of fio the superiority of the choice k = -v/5 is obvious when extrapolating back, i.e. 
to lower scales //. Indeed, this is to be expected from Fig. lb and 2a. 



The resulting ratio Ra{Mz) depends, however, only weekly on A'^' in the interval (200,400) MeV. 
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Figure 2: The same as in Fig. lb, but for fXQ equal to 2 and 10 GeV. 



• If for fiQ = 2 GeV and k = only the b and t quark thresholds (dash-dotted curve in Fig. 
2a) are considered f\, the result is markedly further away from the the exact extrapolation 
than if also the c quark threshold is taken into account. This illustrates the influence of the 
proper treatment of charmed quark threshold even at the energy scales around Mz- 

The influence of the choice of k on the results of the matching procedure has recently been 



investigated, within the approach introduced in [|13[], in Although this matching procedure is 



more sophisticated than that based on the relations (|9|)-(|r|) and is applicable at any order, it 
relates couplants corresponding to different numbers of massless quarks only and thus does not 
concern the difference between the results using the exact expression for the function h{x) and 
those based on its step approximations, as investigated in this note. The conclusion reached there, 
i.e. that as{Mz) obtained by extrapolation from low energy scales changes little (by less than 0.6%) 
when the matching points are varied in some reasonable interval, is therefore not in contradiction 
with the result of this note. 

In summary, the approximate treatment of quark mass thresholds, based on the step-like ap- 
proximation of the QCD /9-function and using the matching procedure of [0] with k = 1, leads 
to an overestimate of a^'^{Mz) extrapolated from low energy quantities, in particular the Rr ra- 
tio, by about 2%. Though not a large effect, it is of the same magnitude as the typical NNLO 
corrections to a^^'^{MS, Mz)- This overestimate of the extrapolated as{Mz) may explain part 
of the discrepancy between the extrapolations based on and deep inelastic scattering. On the 
other hand, the choice of k = \/5, suggested by theoretical considerations, yields extrapolations 
which are in much better agreement with those based on the exact form of the threshold function 
h{x). The analysis of quark mass effects, presented above, holds strictly speaking, at the LO only. 
Nevertheless as both the mass effects the higher order perturbative corrections are small effects, 
it seems reasonable to expect that the main conclusions of the present analysis have more general 

^In these circumstances and for k — the matching point corresponding to charmed quark, given by \/5mc = 3.26 
is below ^0- 
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validity. 

I am gratefull to P. Kolaf and J. Rames for carefull reading of the manuscript and useful comments 
and suggestions. 
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